Role of non-magnetic disorder on the stability of U(l) spin liquid : A renormalization 

group study 
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Recently Hermele et. al claimed that the infrared (IR) fixed point of non-compact QED3 
is stable against instanton excitations in the limit of large flavors of massless Dirac fermions 
cond-mat/0404751 . We investigate an effect of non- magnetic disorder on the deconfined quan- 
tum critical phase dubbed U(l) spin liquid (UlSL) in the context of quantum antiferromagnet. In 
the case of weak disorder the IR fixed point remains stable against the presence of both the instan- 
ton excitations and non-magnetic disorder and thus the UlSL is sustained. In the case of strong 
disorder the IR fixed point becomes unstable against the disorder and the Anderson localization is 
expected to occur. We argue that in this case deconfinement of spinons does not occur since the 
Dirac fermion becomes massive owing to the localization. 
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Recently Hermele et. al pointed out that usual RPA 
treatment of gauge fluctuations pj is not sufficient in or- 
der to examine an instanton effect even in the presence of 
large flavors (N) of massless Dirac fermions in compact 
Qi?Z?3 [3| • They showed that at the infrared (IR) sta- 
ble fixed point of non-compact spinor QED3 in large N 
limit instanton excitations become irrelevant and instan- 
ton fugacity goes to zero. This originates from the fact 
that a magnetic charge has a large value proportional 
to N at the fixed point. The large fixed point value 
of a magnetic charge is due to screening of an electric 
charge byparticle-hole excitations of the massless Dirac 
fermions 3J. The magnetic charge goes to zero in the 
absence of the massless Dirac fermions at low energy Q. 
The larger the magnetic charge, the smaller the proba- 
bility of instanton excitations. As a consequence they 
concluded that deconfinement does exist at least at the 
critical point in the large N limit. In the context of quan- 
tum antiferromagnet stable 17(1) spin liquid (UlSL) is 
obtained^- 

In realistic cases disorder always exists. In the case 
of non-interacting fermions it was shown by scaling ar- 
gument that in 3 spatial dimension the presence of dis- 
order causes a metal- insulator transition But in 1 
or 2 spatial dimension even weak disorder leads elec- 
trons to be localized and only insulating phase is ex- 
pected to exist Q. The presence of long range interaction 
can change the above picture of non-interacting particles. 
Herbut studied role of a random potential resulting from 
non-magnetic disorder on a critical field theory of inter- 
acting bosons via Coulomb interaction 5] . In the study he 
showed that competition between the random potential 
and Coulomb interaction leads to a new charged critical 
point near 3 spatial dimension where a dynamical criti- 
cal exponent z is exactly given by l[5j- In the absence 
of the random potential the charged critical point is not 
expected to exist and only the standard runaway charac- 
teristic is found[5j- In the case of fermi fields Ye investi- 
gated role of disorder on a Chern Simons field theory of 
interacting Dirac fermions via Coulomb interaction in 2 



spatial dimension 6] . In the study he found a line of fixed 
points which is stable in some cases 0. In the absence 
of both the randomness and Chern Simons interaction 
only the runaway characteristic was found as the case of 
bosons. 

In this paper we investigate the role of non-magnetic 
disorder on the deconfined quantum critical phase of 
QED3 in the limit of large flavors of massless Dirac 
fermions. In the concrete we examine the stability of 
the IR fixed point of non-compact QED3 in the presence 
of both the non-magnetic disorder and instanton exci- 
tations. Existence of the stable IR fixed point in the 
absence of disorder is a main difference from previous 
works 0,0- In the previous 

studiesHEI 

there are no sta- 
ble IR fixed points in the absence of disorder as discussed 
above. In the case of weak disorder the IR fixed point is 
found to remain stable and J7(l) spin liquid (UlSL) in 
the context of quantum antiferromagnet is expected to 
survive. The stability against the weak disorder results 
from the existence of the IR fixed point in the absence 
of disorder in 3 space and time dimension. In the case 
of strong disorder we find that it becomes unstable. We 
are led to strong coupling regime where the Anderson 
localization is expected to occur. Owing to the local- 
ization the Dirac fermions become massive. In the case 
of massive Dirac fermions usual RPA treatment may be 
possible. Instantons are expected to be proliferated. We 
argue that deconfinement of spinons is not expected to 
exist in the strong disorder. In addition, we discuss a 
bosonic field theory in the presence of non-magnetic dis- 
order and find a difference in role of disorder on the Dirac 
fermions and bosons respectively. 

First we review deconfinement at the IR fixed point 
of QED3 in the absence of non-magnetic disorder [3 • We 
consider an effective action usually called QED3 in imag- 
inary time 
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Here ij) a is a massless Dirac spinor with a flavor in- 
dex a = 1,...,N and a^, a compact U(l) gauge field. 
x = (r, r) with D — 1 dimensional space r and imaginary 
time r. In the context of U(l) slave boson representa- 
tion of SU(N) quantum antiferromagnet this action can 
be considered as an effective action in the it flux phaseQ. 
In this case the Dirac spinor represents a spinon carrying 
only the spin quantum number 1/2. It is well known that 
QED 3 with non-compact U(l) gauge field has a stable IR 
fixed point in the limit of large flavors of massless Dirac 
fermions 8J. A renormalization group (RG) equation for 
an electric charge is easily obtained to be in one loop 
order 0,11, 



dp 2 



AAe 4 , 
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where A is a positive numerical constant. The first term 
represents a bare scaling dimension of e 2 . In (2 + 1)D e 2 
is relevant in contrast with (3 + 1)D where it is marginal. 
The second term originates from self-energy correction of 
the U(l) gauge field by particle-hole excitations of mass- 
less Dirac fermions. As shown by this RG equation, a 
stable IR fixed point of e* 2 = -At exists in QED 3 . Now 
our question is if the IR fixed point remains stable after 
admitting instanton excitations. Using the electromag- 
netic duality, Hcrmclc ct. al obtained RG equations of a 
magnetic charge g — A and an instanton fugacity y m , 



^L=-(4-D)g-aylg 3 + \N, 



dy n 



dl 
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where a and (3 are positive numerical constants 2] . Ow- 
ing to the last term AA in the first equation a magnetic 
charge has a large fixed point value proportional to A, 
i.e., g* = A A. As a consequence the instanton fugacity 
goes to zero at this IR fixed point. UlSL in terms of 
a spinon (Dirac fermion) and non-compact U(l) gauge 
field is obtained at the quantum critical phase. 

Next we investigate the stability of the UlSL fixed 
point in the presence of non-magnetic disorder. We re- 
consider QED3 in the presence of non-magnetic disorder 
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Here V(x) is a random potential generated by non- 
magnetic disorder. It couples to a spinon density ow- 
ing to the relation of 72] CT =i4 c ff = ^E^=i/]/<t0- 
Here c CT represents an electron with spin a and f a , a 
spinon with spin a. A physically relevant case is that 
the random potential is random only in space but static 
in time. Thus it does not depend on imaginary time, 
i.e., V(x) — V(r). We assume that V(r) is a gaussian 



random potential with < V(r)V(r') >= W5(r — r') and 
< V(r) >= 0H||. Using the standard replica trick to 
average over the gaussian random potential, we obtain an 
effective action in the presence of non-magnetic disorder 
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Here a, a' is replica indices and the limit M — > is to be 
taken at the end. 

Introducing renormalized field variables of ip a = 

2 L Z a ' a^r, we ob- 

e -(5-D- z) l z -l e 2 
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tain renormalized couplings of e 2 = e ( 5 D z ^ l Z a 1 e 2 
and W = e-^- D ~ z ^ l Z- 2 Z w W r . Here z is a dynamical 
critical exponent. Zk, Z a , Z\y are usual renormalization 
constants of a Dirac fermion, gauge field and strength 
of a random potential respectively. A subscript r repre- 
sents "renormalized". Eq. (5) is obtained to be in terms 
of renormalized variables 
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with rescaled space r' 



'r and time r' 



In 



the above we omitted a subscript r for a simple notation. 
Calculating the renormalization constants Zk,Z ai Zw in 
one loop order, we obtain RG equations 
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= (5-z- D)e 2 - AAe 4 , 

= (4 - z - D - X e 2 )W + C(N + c)W 2 (7) 



with positive numerical constants A, \, £, c. Here z is de- 
termined by the condition of e~ 2zl Z^ = e~ 21 Z^ with Z^, 
a renormalization constant of a Dirac fermion in momen- 
tum and Z w , that in energy which gives 



1 + AW 
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with a positive numerical constant A. Our interest is the 
case of D = 2 + 1, i.e., 2 space and 1 time dimension. 
These RG equations basically coincide with those of Rcf. 
[6] in D = 2 + 1 if the term (4 - D)e 2 in the first RG 
equation is neglected. The presence of this term leads to 
the stable IR fixed point as discussed earlier. Precise val- 
ues of the positive numerical constants are not important 
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in our consideration. We solve the above RG equations 
with arbitrary positive numerical constants in order to 
understand a general structure. 

First we check the case of non-interacting Dirac 
fermions in the presence of a random potential. In zero 
charge limit (e — > 0) a RG equation is given by 



^ = {4- z- D)W + ((N + c)W 2 . (9) 
dl 



In D = 3 + 1 (3 spatial dimension) there is an unsta- 



RG equations to the order of i 



ble fixed point of W c = 



with z = 1 



ON+c) 

In D = 2 + 1 of 



to separate a metal and an insulator 
present interest the unstable fixed point becomes zero, 
i.e., W c — 0. Thus only insulating phase is expected to 
exist as discussed in the introduction. 

In the case of interacting Dirac fermions via long range 
"electromagnetic" interaction 9] these RG equations [Eq. 
(7)] show three fixed points in D = 2 + 1; the first is 
W\ c = and e 2 c = with 2 = 1, the second, Wi c = 
1, and the third, W 3c = O(-^r) 



and e 2 ^ = ^ with 



and e\ c = ^ + O(j^) with z = l + O(^r). The first 
is a fixed point of free Dirac fermions which is unsta- 
ble for non-zero charge e 2 0|- The RG flow goes to the 
second, the IR fixed point of non-compact QED3. The 
third is an unstable fixed point. This fixed point does 
not exist in the absence of the gauge interaction. Ex- 
istence of the new unstable fixed point originates from 
the term (4 — D)e 2 of the first RG equation in Eq. (7), 
representing relevance of an electric charge in D = 2 + 1. 
In the case of small strength of the random potential, 
i.e., W < Wsc, the random potential becomes irrelevant 
and the usual IR fixed point (the second) remains stable. 
The stability against the weak disorder results from the 
existence of the stable IR fixed point in the absence of 
disorder in D = 2 + 1, as will be shown below. In the case 
of large strength of the random potential, i.e., W > W 3c , 
the random potential becomes stronger. We are led to 
strong coupling regime where our perturbative calcula- 
tion does not apply. In this case the Anderson localiza- 
tion is expected to occur and thus the Dirac fermions 
become massive. 

In order to see the stability of the fixed points we ex- 
pand the RG equations near each fixed point. Inserting 
e 2 = e\ c + / and W = W 2c + h to Eq. (7), we obtain 
linearized RG equations near the IR fixed point 
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As shown by these RG equations, it is clear that the 
IR fixed point is stable against the weak disorder. The 
stability originates from the finite fixed point value of 
an electric charge, i.e., e| c = t^t. Expanding the RG 
equations near the third fixed point, we obtain linearized 
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dh x 
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dl 
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The second equation shows the instability of this fixed 
point. 

Now we examine an instanton effect on these fixed 
points. Using the electromagnetic duality, we obtain RG 
equations of a magnetic charge g — e~ 2 and an instanton 
fugacity y m in the presence of a random potential 

J = -(5 - z - D)g - aylg 3 + XN, 

^L = U- Z -D- X -)W + ([N + c)W 2 , 
dl g 

^ = (z-l + D- f3g)y m (12) 

with positive numerical constants a, (3 same as those in 
Eq. (3). The term —ay^g 3 in the first RG equation 
is added owing to screening of a magnetic charge by in- 
stanton excitations^- In the case of weak disorder, i.e., 
W < W%c, the IR fixed point of g* — XN remains sta- 
ble as discussed above. Thus the UlSL is sustained by 
the same reason as the case of no disorder. In the case 
of strong disorder, i.e., W > W 3c , the IR fixed point be- 
comes unstable. The strength of disorder becomes larger. 
An important question is whether the magnetic charge 
goes to zero or not as the strength of disorder gets larger. 
At first glance of Eq. (11) the fixed point value of an elec- 
tric charge seems to be sustained. Thus one may conclude 
that UlSL is still expected to occur. In this case, more 
precisely, gapped UlSL may appear owing to the Ander- 
son localization which is expected to occur in the strong 
disorder. But this is an illusion. As discussed above, 
the Dirac fermions are expected to be massive owing to 
the localization. The above calculation cannot apply to 
this strong coupling regime. Further, this phase is not 
expected to be critical any more. In this "insulating" 
phase of spinons screening of internal charge becomes 
negligible. The internal charge is expected to go to in- 
finity following its bare scaling dimension. The magnetic 
charge goes to zero. In this case instanton excitations 
are relevant perturbation and the instanton fugacity gets 
larger to go to infinity. Thus deconfinemcnt of spinons is 
not expected to occur in the strong disorder. 

Now we discuss an effect of non-magnetic disorder on 
a critical field theory of interacting bosons via long range 
electromagnetic interaction. The critical field theory usu- 
ally dubbed scalar QED3 can be considered to describe 
interacting holons via internal gauge interaction in the 
context of U(l) slave boson theory [3, [Tl|. In the scalar 
QED3 a stable IR fixed point called charged or IXY fixed 
point is also found in the limit of large iV^^I as the 
case of the spinor QED3. Here N is the flavor number 
of boson fields. The charged fixed point governs a su- 
perconductor to insulator transition. The charged fixed 
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point is shown to be unstable in the presence of weak 
disorder [Lj- A new stable fixed point is expected to ap- 
pear in association with a random potential. This fixed 
point seems to be related with a Bose glass to supercon- 
ductor transitionjjjl,^!, not a Mott insulator to supercon- 
ductor transition governed by the charged fixed point [l2T|. 
But its nature is not clear. Emergence of the new stable 
fixed point is a main difference between the scalar and 
spinor QED3. This is because even weak disorder is a rel- 
evant perturbation only in a bosonic field theory. The rel- 
evance of weak disorder originates from different scaling 
between bosons and Dirac fermions |LJ|. This different 
scaling is due to a difference in equation of motion, the 
Klein-Gordon equation and Dirac equation respectively. 
At this new stable fixed point the fixed point value of an 
electric charge square is still proportional to inverse of 
the flavor number of bosons, i.e., e 2 c ~ iV _1 ^j|. Thus a 
magnetic charge is proportional to N as the case of the 
charged fixed point. In the limit of large flavor number 
instanton fugacity goes to zero at this new fixed point. 
Deconfinement of boson fields (holons) is expected to oc- 
cur. Bosonic U(l) liquid is sustained at the new fixed 
point associated with non-magnetic disorder. 

In this paper we considered SU(N) quantum antiferro- 
magnet described by N flavors of massless Dirac fermions 
interacting via compact U(l) gauge fields. But a real 
antifcrromagnet has SU(2) symmetry. Thus the flavor 
number of the Dirac fermions is given by N = 2|jj. In 
this case it is not clear whether the present result can be 
applicable. It is known that there exists a critical flavor 
number N c associated with spontaneous chiral symme- 
try breaking (S X SB) in QED 3 R. But a precise value 
of the critical number is still in debate [l^. If the critical 
value is larger than 2, the SySB is expected to occur 
for the physical N — 2 case [3, [l6|. The Dirac fermions 
become massive. In the SxSB phase the massive Dirac 



fermions are confined to form mesons, here, magnons|7j. 
As a result, in the case of N c > N = 2 the U1SL is not 
expected to exist and the present consideration is not ap- 
plied. But there is a cure even in this case. We consider 
hole doping to the U1SL. Doped holes are represented 
by holons carrying only charge degree of freedom^?! • R- e ~ 
cently Senthil and Lee claimed that critical fluctuations 
of holons at a quantum critical point associated with 
a superconducting transition can result in suppression 
of instanton excitations and thus the quantum critical 
point can be described by the U1SL for spin degree of 
freedom ^tJ- This argument is based on the fact that the 
S\SB does not occur owing to critical fluctuations of 
holons. The critical fluctuations increase a flavor num- 
ber of massless fluctuations If a total flavor number 
of massless spinons and holons exceeds the critical value 
N c , the S\SB is not expected to occur [l^. As a 
result the present scenario for the role of non-magnetic 
disorder in the U 1SL has a chance to be applicable. 

To summarize, we showed that the U(l) spin liquid 
is sustained against the presence of weak non-magnetic 
disorder. However, strong disorder leads the fixed point 
to be unstable and the RG flow goes to strong coupling 
regime. In the strong coupling regime our perturbative 
RG does not work any more. Thus more refined cal- 
culation is required. We argued that deconfinement of 
spinons does not occur in this regime since the Dirac 
fermions become massive owing to localization. Lastly, 
we compared the case of Dirac fermions with that of 
bosons. In a bosonic field theory a new stable fixed 
point associated with disorder emerges in contrast with a 
fermionic field theory. At this new stable fixed point an 
electric charge is sufficiently screened and deconfinement 
of bosons is expected to occur. 

K.-S. Kim thanks prof. Han, Jung-Hoon for helpful 
discussions of basic concept in association with disorder. 
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